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Abstract 

The main issue of this work consists in extracting one or several finite 
values for the sum of series involved in perturbation theories. It is sup- 
posed to work for all cases in which two physical parameters are involved, 
and makes thorough use of dimensional arguments concerning these phys- 
ical quantities. Weak and strong coupling expansions are considered as 
the two faces of a formal expression which is the central object of this 
method. This so-called extension procedure is systematic. We apply it 
here to the divergent perturbative expansion of the ground state energy 
of the anharmonic oscillator of quantum mechanics in zero and one di- 
mension, and provide, given a p-order divergent expansion, an analytical 
expression of its estimated sum. This method which is inspired by vari- 
ational procedures provides high degree of accuracy from lower orders of 
perturbation and seems to have remarkable converge properties in a wide 
class of expansions concerning physical observables. 



I. Introduction 



The problem of divergent series has long time been a big concern, especially for 
physicists dealing with perturbative methods. In many cases, the divergence of a 
perturbative expansion arises when the solution of the free part of an equation is 
too "remote" from the exact solution. The choice of this free part is theoretically 
quite open but is actually restrained by the computability of the free solution, 
together with its perturbative corrections. 

When the series converges, the accuracy on the result only depends on our 
ability to compute high order perturbation corrections. Moreover, the conver- 
gence of the series and the fact that the successive approximations tend to the 
right result, might be guaranteed by mathematical theorems. 

When the series diverges, the problem faced by the physicists is one of dealing 
with an apparently worthless result (the truncated series) , despite the fact that 
physical arguments might ensure us that this series actually corresponds to a 
known finite number. The scientific literature is rich of various methods to 
treat this situation: pade approximants [1], variational methods [2,3,4], delta 
expansion [5, 6] with principle of minimal sensitivity [7] or multiple scale analysis 
[8]. 

When facing the situation, one can note the following points: 

1. While dealing with a physical problem, quantities with physical dimen- 
sions are involved, and as a perturbative series expands in powers of a 
dimensionless parameter, there is a first constraint on the way physical 
quantities enter this dimensionless parameter. There is also a second con- 
straint regarding the dimension of the physical quantity expanded itself. 

2. When approximating a physical quantity in powers of the coupling con- 
stant, we sometimes know a strong coupling expansion (in negative powers 
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of the coupling constant) to exist, even if it is not computable. 

These two points are central in the way I propose to deal with a class of 
physical problems involving two physical parameters. The existence of only two 
physical parameters leads to a unique way to build up a dimensionless parameter 
g (say, the coupling constant). 

In the cases where the expression of a physical quantity in term of the physi- 
cal parameters is not exactly solvable, we expect two different expansions of this 
physical quantity in powers of the coupling constant to exist, at least formally: 
the weak coupling expansion (in positive powers of g) and the strong coupling 
expansion (in negative powers of g). Standard perturbative method can give 
us access to the weak coupling expansion up to a certain order p. The case we 
are interested in in this paper is the one when this expansion is divergent. The 
necessary background for our original systematic method to work can now be 
stated the following way: 

the existence of a weak coupling expansion of E, 

oo 

E(g)=J2 E N9 N , (1) 
and of a strong coupling expansion of E, 

oo 

E(g)= g- a/s Y,^g-' lN/s , (2) 

where E is a dimensionless function of the searched physical quantity. The 
parameters 7, a and 5 are crucial ingredients in the correct writing of the srong 
coupling expansion. Their values usually be obtained analysing the equation 
obeyed by the E (see Simon [9] for such a study on the energy of the anharmonic 
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oscillator of quantum mechanics). For our method to work on the calculation 
of E , we thus need to know the following values: 

• the Epf coefficients up to a certain order p. 

• 7, a and S designed to be the smallest possible triplet of integers to sat- 
isfying eq.Q. In the following, this triplet of integer values will be called the 
signature of E(g). 

A good deal of literature has been devoted to variationnal methods using 
principle of minimal sensitivity, which happen to be quite efficient [3]. Some 
other attempts to deal with divergent series through variationnal methods, mak- 
ing clever use of the knowledge of the existence of a strong coupling expansion, 
were successful [10]. The new method we introduce is original in the way it 
considers weak and strong coupling expansions as two faces of a single object 
called an extension. Making use of the formal properties of extensions, one can 
obtain a non-perurbative analytical estimate of E in term of the coupling con- 
stant g, given its signature and a certain amount of coefficients of its divergent 
expansion. Note that the method will also work as an accelerator of convergence 
for convergent series. 

The next section will introduce the extension method for the trivial case 
where the signature is uniform (7 = a = S = 1) so that extensions write as 
expansions in all positive and negative integer powers of the coupling constant. 
The completion of the general method through the actual introduction of the 
notion of signature is performed in the third section. 

The fourth section is devoted to the detailed application of the method to two 
specific cases. It is first experienced on the divergent expansion of the zero 
dimension anharmonic oscillator. The mathematical homogeneity of the inte- 
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grated exponential function allows us to consider the integral itself as depending 
on two parameters with physical dimensions. The extension method will happen 
to reproduce the exact result from the first order of perturbation. As e second 
example, the method is applied to the one dimensional anharmonic oscillator of 
quantum mechanics in the calculation of the ground state energy. It will show 
a very remarkable agreement with the exact result for any value of the coupling 
constant. 



II. Extensions of power series for uniform signa- 
tures 

In the simple case where the signature is uniform (7 = a = 5 = 1), the respective 
weak and strong coupling expansions of E in powers of a coupling constant g 
write: 

00 00 
E(g) =J2 E "9 N and E(g) = g' 1 ]T E N g- N . 

JV=0 N=0 

We will consider these two expansion as parts of a single extension. To do 
so introduce such formal sums as: 

00 

n— — oo 
00 

n— — 00 



00 . 

7 V Ul 7" 

n— — 00 v ' 

The manipulation of negative factorials is made algebraically coherent introduc- 
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ing (—1)! = 0. For all coefficients for which ft do not simplify, one let it go 
to infinity as detailed at the beginning of the third part. When Z is a given 
complex value out of the unit circle, these formal double infinite series have the 
property to be convergent (resp. divergent) in the positive powers of Z (called 
its right part) and divergent (resp. convergent) in the negative powers of Z 
(called its left part), depending on the modulus of Z. 

Note that Z± can be expressed in term of Z using derivatives: 



= 3b. (3) 



They also have the following formal property, 

i>o, (i - zy +1 z t = o , 

from which one can logically expect the left expansion's sum (I) to be in some 
sense "opposite" to the right one (r). The latter expectation can be easily 
checked as long as one performs formal manipulations on power series without 
any consideration about their convergence: Let us define the linear applications 
r and I. 

— oo i 
n=— 1 v ' 

The expression Z'/(l — Z) l+1 , when formally expanded in powers of Z, results 
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in r(Zi), and when expanded in powers of 1/Z, results in — l(Zi). We write: 

r{Zi) ~ 



l(Zi) 



(i - ZY+ 1 
z i 

~ (i _ ZY+ 1 ' 



making sense of our previous remark about riZj) and l{Zi) being 'opposite' to 
each other : r(Zi) ~ — 

In the following, we will combine elements as Zj's to build a vector space. We 
hope any polynomial to have "natural extensions" in this vector space. 

Definitions: 

• Let E p be the C-vector space of all polynomials with a lower than p degree: 

p 

Ae~ p A = J2 A nZ n for A* e C 

n=0 

• Let !B P be the C-vector space generated by extensions Zi for i positive and 
lower than p: 



V e S p =^V = J2V 1 Z.for V 1 G 



II. 1. Extensions of a polynomial 

An element of S p will be considered as an extension of a p-order polynomial if 
its restriction to zero to p powers of Z identify to the given polynomial. We 
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show now that this properly defined extension of a polynomial is unique. 
Definition: 

Let A and V be respective vectors of S p and S p , then V, once expanded as: 

oo 

v= Yl v ^ zn > 

n= — oo 

is said to be a p-order extension of A if: 

V n = An for < n < p . 
Moreover, if V = Y^=o ^i, from the definition of the Zi's, we have: 



^— ' i!(n — in 

i=0 y ' 



We will now proceed to show that such a p-extension of a given polynomial is 
unique through the introduction of a specific linear mapping between S p and S p . 

Theorem 1: 

For given p, let ip and ip* be the following linear mappings: 



if: 3" — E p 



EP i\ t 1\n—iy 

i=n n\(i-n)\ V ^ Zv « 



• 77 > 77P 



EP n! 7n 

n— i i!(n— i)! 
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Then ip and ip* are inverse linear mappings. 



Proof: 



i=n n\(i—n)\ 



(-i)»-V*(^i) 



EP y^P 



fc=i n!(»-n)!(fc-*)! 



(-l) n - l Z k 



Using the following identity for k < p, where S is the Kroeneker symbol: 



>From the existence of the reversible mapping ip between and S p and 
the fact that {^™}o<n<p is a base of S p , it comes that {i?i}o<j<p is a base of H p . 

Theorem 2: 

Let i € 5 P , then € 5 p is the only extension of A. 

Proof: 

Let us first recall that: 



such that by construction, for < i < p, Z% is a p-order extension of ip*{Zi). 




we prove that for < n < p and < i < p: 



tp*(tp(Z n )) = Z n and ip{p{Zi)) = Z, 
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If now A = J2n=o A n Z n , one can write <p(A) — Y^!=o ^ where: 

A * = E ,r tl ■ ( 4 ) 

' n (z — n ! 

n=0 v ' 

>From the reversibility of if, one has if (if* (A)) = A, and we can establish that: 

^) = ELo Al z i ■ 

As Zi is an extension of <p*(Zi), it follows that if (A) is an extension of A for 
any A S S p . The reversibility of y> guarantees the unicity of this extension. 

Example: 

For p = 2, let us find the extension of the 2-degree polynomial: 

A= 1 + 2Z + 3Z 2 . 

The coordinates of its second order extension are given by the relation Q : 

A = A =1 
A 1 = -A + A 1 =1 
A 2 = A -2A 1 + A 2 = , 

so that: 

• if {A) — A Zq + A 1 Z\ is the extension of A; it also writes: 

oo 

<p(A) = ( n + 1 ) zn > 
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and is actually a second order extension of A, as expected. 
• r(ip(A)) = r{Z Q ) + r(Zi) ~ 1/(1 -Z) + Z/(l - Zf = 1/(1 - Z) 2 is what one 
could reasonably expect as a sum for a divergent series behaving like A does for 
lower orders. 

Remark: 

This example illustrates an important case where the extension of a p-order 
polynomial (in S p ) actually belongs to S p_1 , that is A p = 

The latter remark will be of great use in the following introduction of a 
parameter in our method. It will also be verified as a test for the reliability of 
the method in the example of application given in the fourth section. 

II. 2. The introduction of a parameter: uo 

We will show that the simple introduction of a complex parameter u> allows us 
to build as many S p vector spaces which means as many p-order extensions of 
a given polynomial as non vanishing values for u. There is then no reason why 
restraining to the previous u — 1 case. At the end, the problem will just remain 
to find the most appropriate value of u>. 

>From now, for a fixed u, E p and S p will be redefined: 

p 

AeZ p A = A n {uZ) n ioruj G C* and A n G C 

P 

v e s p =>v = J2 yi i uZ )i for w e c * andyt e c 

i=0 
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where (u>Z) n = u> n Z n , and (uiZ)i is defined as follows, 
p now acts on vectors as (w2) n : 

Theorems 1 and 2 still hold as long as to does not vanish and the method now 
produces an extension <p{A) of A, depending on the value of u>: 

A = ±A n Z" = ±^(uZr , 

n=0 n=0 

*w = E ^ p((<^n = E A » (^)* . 

n=0 i=0 



where 

y ' ^ n\(i - n )r ' u) n 

n=0 y ' 



When trying to extend a given p-degree polynomial A, one now has to deal 
with as many extensions as possible values for u>. The last remark of the preced- 
ing section suggests that we could use this freedom on the value of the parameter 
ui to make A p (tu) vanish so that we could just keep p-order extensions of A that 
lie in S p _i. A p (uj) being a p-degree polynomial in u>, one expects u> to take at 
most p values (A p {uj) = 0) so that a vector of S p will have at most p different 
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extensions. If vo is a root of A p (lu), it solves: 



■777^ 



n!(p — n)! 



(5) 



and the associated extension writes: 

p-i 

<p(A) = 2 AV) (n7Z) 4 . (6) 

i=0 

The problem now remains to choose from all these extensions. 
II. 3. A criterion for the choice: u> = w 

The aim of this work is, given a divergent series up to order p, to obtain an 
estimation of its complete sum. Let us start with the method applied to the 
(p — l)-order. In order to establish a convenient criterion for the choice of the 
"better" ui between the roots of we will look for the one which makes 

the better prediction for the forthcoming term of the series: A p . 
To do so, let us suppose that A p_1 (w) and A p (u>) have a common root w, the 
following properties come: 

• From the equation ©, the tn-associated (p— l)-order and p-order extensions 
identify. As a consequence, the forthcoming term predicted by the (p— l)-order 
extension is the exact value A p . 

• From the definition © of A p (lu), one can establish the following relation: 



so that if w is a common root of A p 1 (oj) and A p (lu), it has to be a double root 
of Ap(lu). 



= —pur 
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The goal is to achieve the best estimate of the complete sum of the series 
through its extension at a finite order. It seems natural to be interested in 
the case where two successive extensions (orders p — 1 and p) are equal, thus 
at order p, double roots are our best candidate for the choice of w. In the 
absence of double roots, we decide to look for the closest couple of roots - in 
terms of distance in the complex plane - and select one of these. In doing so, 
we expect our selected root to be also close to a root of the previous estimate. 
As a consequence, these two successive extensions are hoped to be close to each 
other. The applications of the method in section IV will corroborate, through 
numerical observations, these general remarks. 

III. Signature of series in physics 

For the moment, to keep reasonable, one has to expect most of the series en- 
countered in physics not to fit in this frame in the sense that the predictions for 
the forthcoming terms of a development are far from the right one (and that 
the roots are far from each other) . Let us consider for example such a "misfit" 
expansion: 



which expected large order behaviour can obviously not be found as a combina- 
tion of the Zi's as defined. The fact that forthcoming terms of such a series will 
be very badly predicted by our method can also be pointed out when noticing 
that an extension of this truncated series has to vanish for even and negative 
values of n (as factorials of negative integers diverge). The Zi's as defined can 
obviously not satisfy such a condition. This section will take into account this 
kind of difficulty (now to be considered as a clue) when introducing the notion 
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of signature relying on the knowledge we have of the form of the strong coupling 
expansion. 

Given a finite series of complex numbers, we now have at our disposal a 
method which permits to extend it infinitely together in positive and negative 
powers of the coupling constant Z . Moreover, the right and left parts of this 
extension can be summed formally. Some of the series we are interested in are 
characterized by cancellations of an infinite set of coefficients (e.g. A n = for 
n even and negative. We shall later concentrate on how the occurrence of the 
vanishing coefficients is in related to the form of the strong coupling expansion. 
Since no combination of the Zi's can reproduce this kind of behaviour, we ask 
ourselves now what to do then when we know the expected series to have such 
a typical property? In this paper, we shall consider such cases where the occur- 
rence of vanishing terms is regular. Let us first consider the following double 
infinite series where by construction all even and negative powers of the coupling 
constant do vanish: 

T. t T.Z- where T. - tl Z "J^Z ■« (7) 

n= — co ^ ' 

As previously mentioned, it will be of great use to reconsider the definitions of 
the factorials in order to make sense with expressions as (— 2)!/(— 1)! = — 1: 

Definition: 

• (n — l)!n = n! for n € Q*, 

• (— 1)! = fl is infinitely large, 

• n! = T(n + 1) for n e [0, 1] n Q. 

The following properties come: 
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• n\ (—1 — n)! sin(n7r) = tt for n £ Q and n ^ Z, 

• n!(-l-n)!(-l) n = for n e N. 

Using these properties, the only non- vanishing terms of T are: 

T]v = ffl^i)Z2!(- ir and T-^-S&i-ir for » eN 

When looking at the left part of T, these cancellations are typical of what has 
been called the signature of T, as it points out all the vanishing terms of the 
expansion it corresponds to. 

Definitions: 

• If T is such that its non- vanishing terms are Tsn and T_ a _ 7 jy for TV £ N, then 
it will be said to have a j\a\S -signature. A typical example of this signature is 
provided by: 

_ (-1 -f)!(-l + g±a)I 

n 

For a given signature, this T will be called its associated transformation. 

As a matter of example, the above series (0 has a 2|l|l-signature and the 1|1|1- 

signature considered in the preceding section was simply called uniform. 

• Let us now define the multiplication T * A: 

oo oo 

A= A " ZU T*A= T n A n Z n 

n—~oo n=— oo 

then, if A has a uniform signature, T * A will have a j\a\S -signature through 
the action of the transformation T. 
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• If B is an extension of A, then T * B is a T -extension of T * A. 
III.l. The case of non-uniform signatures 

As it is formatted to give rise to extensions with uniform signatures, the method 
developed in the preceding section obviously needs a few adjustment to fit the 
cases where the signature is known not to be uniform. Suppose we wish to 
extend a polynomial B known to be associated to a non-uniform signature. The 
method now consists in "removing" its signature multiplying it by T -1 , and let 
it be extended (through the action of ip) into an extension with uniform signa- 
ture. Then, we just have to put back the proper signature multiplying the result 
by the associated transformation T. We thus make sure that the T-extension 
obtained has the expected signature (respecting the expected vanishing terms) , 
and is actually a T-extension of the original polynomial. 

If B is associated to a 7|a|<5-signature, then it expands in powers of Z s in 
the weak coupling limit and writes up to order p: 

v 

B = J2 B Sn Z 5n 

n=0 

Where the connection with equation (TIJ is obtained replacing g with Z 5 . B can 
then be rewritten T * (T^ 1 * B) where T~ l * B can be associated to a uniform 
signature: 

p 

T- 1 *B = Y,T Sn 1 B Sn Z 5n 

n=0 

The extension of T _1 * B can now be found using the method of the second 
section. </?(T _1 * B), its extension, will have a uniform signature. The latter 
product with T ensures us to recover the proper 7|a|<5-signature for the T- 
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extension of B: 



T*ip{T~ 1 *B) is a T — extension of B 



Definition: 

In order to achieve the previous manipulation ip and ip* have to be formally 
updated in replacing Z by Z 5 in their definitions so that they remain inverse 
mappings. The formal expressions then arising are made explicit as follows: 



Using this procedure, we can now obtain the unique extension of any polynomial 
given in addition to its coefficients, its non-uniform signature. 

Our goal is at the end to obtain a non-perturbative expression corresponding 
to the physical quantity corresponding to B up to a finite order. To do so, we 
just have to calculate the sum whether of the right part or of the left part of an 
extension with a non-uniform signature. Both these sums will reveal to equate 
up to a simple multiplicative term. If one of them is convergent, the other is 
not and we will have to consider this second one as formal. 

III. 2. Summation 

In the general case of a 7|a|(5-signature, left and right parts of a T-extension B 
of B write: 




oo oo 



r(B) = B SN Z m and 1(B) = ^ B^ a ^ N Z 



— a — 'yN 
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The p-order T-extension B writes T* A where T is the transformation associated 
to its signature, and A = </?(T _1 *B) is a linear combination of vectors ((wZ)*),: 



_ (-1 -f)!(-l + 



p-i 

A = ^ ((wZ) s ), where w satisfies A p (u;) = 

i=0 



Now recall from equation <EJ| a useful relation between extensions: 



{(uZ)% = ±Z 5i {^y{{uZ)% (8) 



This property suggests to rewrite T * A as: 

T*A = Y t A%)±Z*{^ g ) i [T*((cZ) 5 ) ] (9) 

Then, from the sums of T * ((ojZ) s )o, the sums of T * A will follow for any 
extension A. But as (Z s )o = Zq and T * Zq = T, we simply have to sum the 
left and right parts of T. 

The summation of T : 

l l J J — 2^N=0 1 ~ a ~~fN ^ — 1^N=0\ L ) jvl Z 

r(T) = ESU^Z™ = EN=o(- l ) N i ~ 1+(a m N]h)] z6N 

Let us now introduce integral expressions for fractional factorials where N and 
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v are respectively positive and non zero positive integers: 



{N/v)\ = v 



dtf 




Jo 

If we accept to permute the sum and the integral without wondering whether 
they diverge or not, left and right parts of T can be formally summed: 



The link between left and right parts of T comes from a simple change of variable 
in the integral: 



The analytical properties of I and r as functions of Z are the one of multi- valued 
functions. It has also to be noted that, depending on the values of 7, a and 5, 
the integral expression does not always make sense as conditions on its existence 
might not be satisfied. Nevertheless, as 7 and S are positive integers, only a 
strictly negative value for a can make the integral diverge around t = 0. If this 
situation occurs, then the first terms of the expansion l(T) and r(T) (in respec- 
tive powers of Z~ J and Z s ) will not be incorporated in the integral expression. 
As only a finite number of terms are concerned, this will be a complication but 
not an obstacle to the summing of the series. The application of the method 
given in the fourth section will illustrate this situation. 

The summation of T * A 




and 




7/(T) ~ Sr(T). 



(10) 
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The summed expression for r(T * A) is now from ®: 



(11) 

- 7 / 00 fe-( t /^ T - t5 (^) Q Ef=o^M^(-^ 

where one has successively performed all the derivatives and then a change of 
variable. The first expression has, by construction, the correct expansion in 
powers of Z s up to order p. The relation 1)1011 between left and right parts of 
T * A still holds as I and r are linear mappings: 

yl(T * A) ~ Sr(T * A) 

The formal manipulation performed to obtain the second expression of l(TT|l is 
such that its expansions around Z ~ oo and Z = will respectively identify to 
^l(T*A) and r(T* A). This last point is crucial in the way it it guarantees that 
the form of the left and right parts (the strong and weak coupling expansions) 
will be properly recovered. 

The appendix will concatenate the step by step expansion procedure to get an 
extension and a sum from the knowledge of a p-order series and its signature. 
The next section will illustrate the method on the divergent perturbative series 
of quantum mechanics in zero and one dimension. 

IV. Applications of the extension method 

The extension method will now be applied to the case of the zero dimensional 
anharmonic oscillator, which happens to be a simple integral. It will then be 
tested on the calculation of the ground state energy of the quantum mechanics 
anharmonic oscillator. These applications will follow the extension procedure as 
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it is concatenated in the appendix (the signature, the extension and the sum). 

IV. 1. The zero dimension anharmonic oscillator 

We consider the following integral: 

/>OC 

E(m,X) = V2 dxe- mx2 ' 2 ~ Xxi /\ 
Jo 

which, when expanded in powers of X/m 2 gives rise to a divergent series. 
The Signature: 

A dimensional analysis on the base of the zero-degree homogeneity of the expo- 
nential function gives: 

[E] = [x] = [m- 1 ' 2 ] = [A- 1 / 4 ], 

where m and A are to be considered as the two parameters with a physical 
dimension. Then E will expand both in powers of X/m 2 in the weak coupling 
limit, and in powers of (A/m 2 ) -1 / 4 in the strong coupling limit: 

oo , 

E(X)=m-^^E N (±) N , 

OO 

E(X) = X-y^E N ( ™ )». 
Note that -E'(A) is divergent and E(X) is convergent. 

The signature in then given by (7 = 2, a = 1, S = 4). Taking X/m 2 = Z and 
m = 1, right and left expansions of such an expression respectively rewrite as: 

00 00 

N=0 N=0 
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exhibiting a 2|l|4-signature. 



Given the first coefficients of the right expansion, let us apply the extension 
procedure to the first order of this divergent series (now called B) : 

B = B + B 4 Z 4 + ... 



where 



and 



/>oo 

B Q = V2 dx e - mx2/2 = 
Jo 



B 4 = -V2 dxe~ mx I 2 —=--^. 

J 

We now need the transformation associated to the 2|l|4-signature: 

(-1 -f)!(-l + l±g)! 
SI 

such that: 

r- 3 r- 

T = V¥ and T 4 = ~^Vk. 

The first two coefficients of the uniform expansion (freed from its signature) to 
be extended write: 

A) = T^Bq = 1 and A 4 = T^B 4 = 1. 



The Extension: 

The parameter uu has to satisfy the equation: 

A\uj) = -1 + — A <^=> tu 4 = l 
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We have now all the necessary material to give an estimation of the whole sum, 
given its first two terms and its signature. 

The Sum: 

The first order estimation of the integral through the extension procedure is 
then: 

/■CO 1 

2/ At L e -w?-t\ 

Jo Z 

which also writes (taking t = y/mZ xj\/2 and Z = (A/m 2 ) 1 / 4 ): 

/>oo 

V2 dxe- mx2 / 2 - Xxl /\ 
Jo 

This an exact result. 

As a matter of comparison: in the case of the calculation of this integral, 
the usual variationnal method with principal of minimal sensitivity happens to 
provide a non-exact finite result from any of its p-order divergent expansion 
[2,4] while the extension method gives the exact result from the first order of 
perturbation. It is now straightforward to notice that any integral of the form: 

/•OO 

/ dxx^e-*"- 2 '*' 
Jo 

will also be exactly estimated form its first order of expansion in powers of Z s 
using the extension procedure. This class of integral can thus be considered as 
canonical regarding the extension method. 

IV. 2. The ground state energy of the anharmonic oscillator 

The development of the vacuum energy E of the anharmonic oscillator of quan- 
tum mechanics is known to be divergent when expanded in powers of the cou- 
pling constant A. Many papers tried with success to deal with this divergent 



24 



series, sometimes with high level of accuracy and/or complexity [3, 7, 10, 11]. 
The Lagrangian of the model writes: 

L(m 2 ,A) = \{d t <5>f -\m 2 & _ 1 A$4 

A dimensional analysis of physical quantities gives: 

[L] = 1 ; [t] = -1 ; [$] = -1/2 ; 

[m] = 1 ; [A] = 3 ; [E] = 1 . 

The vacuum energy for the free Lagrangian E(m 2 ,X ~ 0) and its perturbative 
corrections in powers of A are computable, but one can also imagine another 
perturbative expansion in powers of m 2 around the free Lagrangian L(m 2 ~ 
0, A). This latter strong coupling expansion, although proved to exist and be 
convergent [9], is not computable with the usual technique of perturbation as 
one cannot solve the free equation. Nevertheless, from dimensional analysis, we 
know that these two expansions respectively write: 

oo , 

E(X) =mJ2E N (-3 ) N 

N=0 m 

for the divergent expansion around A = and, 

oo 2 

E(X) = AV3 J2 e n ( - )- 

JV=0 

for the convergent expansion around m 2 = (or A ~ oo). 

Without loosing any dimensional information, one can get rid of fractional pow- 
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ers writing X/m 3 = Z 3 and m = 1, then: 

oo oo 

E = E N Z 3N and E = Z Y E N Z~ 2N 

N=0 N=0 

Considering both these series as right (E) and left (E) parts of the same object 
we can recognize an extension with 2| — l|3-signature. 

We can now use the extension procedure to sum the known divergent expansion 
of E to a finite order p on the grounds of a 2| — l|3-signature (7 = 2, a = — 1,6 = 
3). The perturbative series of the ground state energy E to be extended will be 
called B. A systematic procedure computes the coefficients of the ground state 
energy (B 3n ) [12]: 

1 3 , 21 „ 333 „ 30885 12 916731 

B = — I Z 3 Z 6 H Z 9 Z 12 H Z l0 + . . . 

2 16 128 1024 32768 262144 

Left and right parts of its extensions are to be respective approximations of \E 
and E. 

IV.2.1. The choice of the "best" root 

The search for a p-order T-extension goes through the choice between p roots 
of A p {tu): (Ji, uj 2 , ■ ■ , u p . The discussion of the section II. 3 suggests to consider 
the root which is as close as possible to another one of the same order (a similar 
comparison with the preceding order roots reveals to be equivalent). In order 
to measure this property of a given root un to be close to one or several other 
roots of the same order, we introduce its weight: 

weighty) = V — ^ 
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The following array will give the roots of A p (uj) for p = 2, 3, 4, and 5 together 
with their weights (in brackets): 



p cut (weight) 

2 1.062 (191) 1.135(191) 

3 1.053 ± 0.06H (113) 1.189(90) 

4 1.043 ±0.102i (135) 1.071 (216) 1.236 (79) 

5 1.034 ±0.14i (135) 1.073 ±0.033i (355) 1.279 (71) 

For every order, the root associated to the biggest weight are underlined (the 
choice for the underlined root of the second order is arbitrary as both of them 
have the same weight). 

As A p (cu) is in fact a polynomial in powers of 1/w 3 , it has p triplets of roots, is 
why we decided to consider only one root in each triplet: the one lying around 
the real axis. This choice will not affect the right part of the extension. 

In order to test our selected roots, a similar array will show the related 
estimations for the forthcoming term of the series (as the p-order T-extension 
provides an estimation for the known value of B^ p+ ^). The estimations associ- 
ated to the chosen roots will be underlined and the last column will show the 
exact value: 

P B 3p+3 (est.) B 3p+3 (exact) 

2 0.3223 0.3211 0.3252 

3 -0.9411 ±0.0035 -0.9552 -0.9425 

4 3.5034 ±0.008 3.4948 3.4269 3.4970 

5 -15.6620 ± 0.0049 -15.6165 ± 0.0034 -16.2018 -15.6208 

The relative precision on the results in the complex plane for the chosen roots 
for respective orders p = 2 to 5 are 9. 1(T 3 , 4. 10" 3 , 6. !CT 4 and3. 1CT 4 
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As far as I know, the criterion governing the choice for the best root still works 
systematically for higher orders of perturbations. 

IV. 2. 2. The vacuum energy in the strong coupling limit 

The vacuum energy of the massless theory identifies to its strong coupling limit: 
A 1 / 3 ^ where Eq is approximated by %B\. One can verify that the choice for 
the best root made earlier still holds. The following array shows the estimated 
values and their relative precisions compared to a known value [11] of Eg for 
order p = 1 to 5: 



p 


\B\ 




1 


0.420139 


(2.10- 3 ) 


2 


0.4205216 


(1.10" 3 ) 


3 


0.4208109 ± 0.0000952 


(2.10- 4 ) 


4 


0.4207976 


(2. 10~ 5 ) 


5 


0.4208087 ±0.0000033i 


(1.10" 5 ) 


E 


0.4208049 


(exact) 



The accuracy obtained for the first order (where only the two first terms of B are 
involved) is remarkably high (0.2%) and increases with the order of perturbation. 
We could also check that the choices for the best roots systematically hold for 
all the coefficients B n , and also for the sum itself. 

IV. 2. 3. The summation for finite coupling 

We also dispose of an analytical expression of the vacuum energy in the case of 
a finite coupling (for any value of A = Z 3 and m = 1). In order to compute 
numerically this integral, the zeroth order term in powers of Z 3 of the integrant 
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of the right sum (which we know to correspond to 1/2) will have to be extracted 
from the integral as it is not summable around t = 0. As an example, we give 
the analytical expression of the first order estimation for E(X): 

where w ~ 1.09954. 

The following arrays show an estimation of E(X) for A = 2, 4 and 8 to be 
compared with the exact value [10,11]. The corresponding relative precisions 
are in brackets. 



V 


A = 2 




1 


0.6957043 


(7.10- 4 ) 


2 


0.6950365 


(2.10- 4 ) 


3 


0.6951801 ±0.0000363i 


(6. 10" 5 ) 


4 


0.6961732 


(3.10- 6 ) 


5 


0.6961768 ±0.0000008i 


(1.10" 6 ) 


E{2) 


0.6961758 


(exact) 
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p 


A = 4 




1 


0.8030005 


(1.10- 3 ) 


2 


8035964 


C\ 10 _4s l 


Q 
O 


u.ouo i * / ox u.uuuuoyoi 


(°- 1U J 


1 


U.oUo i 000 


(r in- 





n sn^7797 -i- n nnnnm cw 
u.ouo iiZi ± u.uuuuuiyz 


I 4 - iU ) 


£(4) 


0.8037707 


(exact) 


P 


A = 8 




1 


0.9504142 


(1.10- 3 ) 


2 


0.9501856 


(4.10- 4 ) 


3 


0.9505784 ±0.0001160i 


(1.10- 4 ) 


4 


0.9515597 


(1.10" 5 ) 


5 


0.9515723 ±0.0000034i 


(5.10- 6 ) 


E(8) 


0.9515685 


(exact) 



As expected, the selected roots provide systematically the best approximations, 
which reveal to improve greatly with the order of perturbation. 
One can also test successfully the procedure on other observables of the model, 
such as the mean value of < $ 2 > or the excited states of energy. The same 
relative precisions is achieved in these calculations. 

Eventually, one can reasonably expect the method to work as well for any 
x 2N potential in one dimension quantum mechanics. 
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V. Conclusion 



This paper does not provide any proof of the convergence of the extension 
method despite the logic inherent to its construction. Nevertheless, in the case 
of the anharmonic oscillator, it reveals to have a highly predictive power for all 
positive values of the coupling constant. It also makes a crucial and automatic 
link between weak and strong couplings and provides an analytical expression 
for the sum of some divergent perturbative series. 

Regarding the problems of divergences arising from quantum mechanics and 
more generally from domains of physics facing divergent expansions involving 
two physical parameters, the method exposed here, besides its simplicity, has 
the advantage to be systematic. It is also hoped to provide an excellent approx- 
imation of a physical quantity (for any value of the coupling constant) from the 
first perturbative correction. It opens perspectives for field theories where only 
lower orders of perturbative corrections (expected to be divergent) are available. 

The power of the method obviously comes from our ability to take into 
account physical information about dimensional quantities. If we want to make 
it efficient for more complex models, or prove that other classes of series suits the 
method, we will have to extend it to models involving more than two physical 
parameters. To do so, we will manage to exhibit, in the same spirit as in this 
work, constraints arising from dimensional information. This will be the subject 
of a second part of this work to be published. 
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Appendix: the extension procedure 

Given respective left an right forms of expansions: 



J2 B- a - T a Z- a -^ and J2 B 5n Z s 



Z5n 

n=0 n=0 



We now summarize the method developed in this paper to estimate the value 
to be associated to the known divergent truncated series: 



]T B Sn Z Sn . 



n=0 



The Signature: 

We are in the general case of a 7|a|<5-signature, the associated transformations 
is: 

(_l-2)!(-i + s±2)! 



The Expansion: 

A Sn =T 5n 1 B Sn 0<n<p 

AP{u) = => uj = w 



B n =T n A n ne Z 
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The Sum: 



00 r oo 1, P- 1 1 



which can also be written: 



The Weight: 

Between {wi,o;2, u p } roots of ^4 p (w), w is the one with the biggest wei 



weight^) = V - — - 
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